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Abstract 

We calculate thermodynamic properties of a disordered model insulator, 
starting from the ideal simple-cubic lattice {g = 0) and increasing the disorder 
parameter (7 to ^ 1/2. As in earlier Einstein- and Debye- approximations, 
there is a phase transition at Qc = 1/2. For g < the low-T heat-capacity 
C ^ whereas for g > gc, C ^ T. The van Hove singularities disappear 
at any finite g. For g > 1/2 we discover novel fixed points in the self-energy 
and spectral density of this model glass. 
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The anomalous thermal properties observed in dielectric glasses [0, 0, has 
prompted numerous theoretical and experimental studies on the subject 0]. 
Arguably, the best well-known model consists of the phenomenological two- 
level states (TLS) proposed by Anderson[^] some time ago. Although the 
physical nature of the TLS remains elusive the theory does predict linear 
specific heat at low-T and enhanced ultrasonic attenuation, both universal 
properties of glassy materials [^. An alternative microscopic theory entirely 
based on phonons was subsequently advanced by one of the authors and his 
collaborators 0. It relies for its results on the interplay between disorder, 
parametrized by a dimensionless and anharmonicity, parametrized by a 
dimensionless ^ (the Gruneisen parameter) that is ultimately taken to 0. The 
solutions of two "toy" models, the one based on the unperturbed Einstein 
model of the phonon spectrum in the "crystalline" phase and the other, on 
the Debye model 0], both exhibited the following features: at a critical 
threshold Qc the system transitions from a "disordered crystal" phase with 
partial long-range order (LRO) to a "glassy" phase, characterized by a linear 
specific heat and a divergent Debye- Waller exponent-i.e. by a total loss of 
LRO0. 

Here, for the first time, we have extended this type of model disorder cum 
anharmonicity into the study of a bone-fide crystal and have obtained closed- 
form solutions in the simple-cubic lattice. This approach is the most realistic 
so far, given that initially the phonon density-of-states of the ideal crystal 
exhibits the van Hove singularities (vHS) thought to be characteristic of a 
crystal with LRO. The results are also much richer than in the toy models 
and can, in principle, be contrasted or compared with molecular dynamics 
performed on the identical simple-cubic lattice. 

Principally our findings are as follows: In the disordered crystalline phase, 
where the vHS might be thought to weaken with increasing g and ultimately 
disappear, instead they soften and disappear immediately at any finite g > Q-, 
despite the persistence of LRO and of a finite Debye- Waller exponent up until 
9 > 9c- There is an exothermic first-order phase transition as g is increased 
above gc = 0.5. Nevertheless the elastic properties in the glassy phase are 
intimately related to those of the crystal, as we shall show. Additionally, in 
the glassy phase g > gc, we discovered that the spectral density function p 
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exhibits fixed points: frequencies at which the spectral density is independent 
of g. This feature bring us close to the ultimate goal of being able to express 
a microscopic law of corresponding states for glassy materials. The specific 
heat exhibits a Debye law throughout the crystalline phase, up to gc- The 
exact coefficient B{g) is given in closed form below. Beyond gc the low-T 
specific heat is linear in T, with coefficient A{T) also given in closed form 
below. 

THE HAMILTONIAN. In an effort to make this paper self-contained 
we review here the main features of the model. The starting point is the 
Hamiltonian of the reference crystal 

Ho = Y.^k(aia, + ^^ {h = l), (1) 

where k contains the wavevector and polarization of the phonons. The opera- 
tor al (ofc) are the familiar creation (destruction) phonon operators. Disorder 
is introduced in the model via the rather general phonon scattering term 

Hi = jY.M(^^^')V^^QkQk', (2) 

k,k' 

where Qk = alfc + ctfc and the M{k, k') are random scattering amplitudes 
which depend on the location and type of disorder. Anharmonic effects are 
introduced next in order to stabilize the spectrum, since Hq + Hi possess an 
unstability threshold where the phonons become over coupled [^, ^. A simple 
quartic term suffices for this purpose, 

'^'^k^k'QkQ-kQk'Q-k' 1 (3) 

where ^ is the Gruneisen parameter, typically a small quantity, expressed 
here in dimensionless form. N is the number of normal modes and ujq is a 
characteristic frequency of the reference crystal. We model M(A;, k') by a set 
of random phases (g/^/N) exp[6'(fc, k')], where the {6{k,k')} are random and 
distributed uniformly in [0, 27r] and g measures the strength of the disorder 
{g"^ is proportional to the concentration of defects in the crystal). Then, 
after a symmetry-breaking transformation is performed and the limit ^ —>■ 
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is taken, we find the ground state energy is a discontinuous function of g at 
Qc = 1/2. Tlie (i?/nam«ca/ Hamiltonian|^, |^ governing lattice vibrations is, 
liowever continuous, and independent of ^ in the hmit. It is: 

H = "^ujk (alak + ^) + 7 XI V^ki^k'M{k, k')QkQk' 



k,k' 



+0 U-oPU-o E^fc^fcQ 



(4) 



where is the Heaviside step function. The model is characterized by the 
existence of a 'critical' disorder concentration Qc = 1/2. In the disordered 
crystal phase {g < 1/2), the anharmonic corrections can be neglected, while 
in the "glassy" phase {g > 1/2) they give rise to an additional quadratic 
contribution responsible for an anomalous shift in the density-of-states to 
low frequencies. This, in turn, results in C ~ T rather than at low 
T. This remarkable mechanism is quite general, transcending the particular 
dispersion uj{k) of the reference crystal. 

The free energy of the model can be expressed as an integral over the 
coupling constant. 



duj[ po{uj) + Pi{lj) + P2{lj) ] log(2 sinh(/3u;/2)). 



where po{uj) is the density of states of the reference crystal. 



Po{uj) = — u Im 

71 



FBZ (27r)3 u;^ — ujj. 



(3 = l/ksT 
(5) 

(6) 



and the spectral densities pi{uj) and p2(^^) are the contributions from disorder 
and anharmonicity0: 



Pi(^) 



TT g- 
^{g - 1/2) 
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/{2g)-i A'^ duo 

and R{uj, A), /(tu. A) are the real and imaginary parts of the the self energy 
Z{uj, A) which obeys the trascendental equation^ 



-I{uj,\)R{u,\)] (7) 
-{2\g-l) A[_/(^,A)], (8) 



7( l\ (^^)V 



cu2-cui(l + + Z(cj,A)) 



(9) 



with = {2\g - 1) ^{\g - 1/2). 



4 



RESULTS FOR THE SIMPLE-CUBIC LATTICE. The integral 
(ID can be performed analytically if the dispersion relation are ul = ti;o(3 — 
cos(A;i) — cos(A;2) — cos(A;3)). After inserting this dispersion into Eq.(|^) and 
taking the continuum limit, the equation for the self energy Z{uj^ A) can be 
written implicitly as, 

W [3, t{u)] = (^1 - 3^(1 + + Z{u))^ + + + z{lo))^ , 

(10) 

where the argument t{uj) is given by, 

. . 3(l + + Z(^)) 

^ ^ 3(l + + Z(^))-(^/^o)' ^ ' 

and Vr[3,z] is the generalized Watson integral evaluated by JoyceQ, who 
expressed it compactly in terms of complete elliptic integrals of the first 
kind§: 

iy[3, z] = l—^K{k+)K{k.) , (12) 
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k% = (1/2) ± (l/4)x2V4 -X2- (l/4)(2 - X2) VI - ^2- Figures 1 and 2 show 
the real and imaginary parts of Z{uj,g). The imaginary part I{uj,g) shows 
the not unexpected systematic broadening of the bandwidth with increasing 
disorder. R{uj, g) exhibits a change in curvature from convex to concave at 
uj = a.s g is increased from below to above gc = 0.5. Correspondingly, the 
slope of I{uj,g) vanishes at = for all g < 1/2 but acquires a finite value 
for ^ > 1/2. 

Z{uj,g) can be expanded at both the low and high frequencies. It is also 
seen to possess a low-frequency fixed point in the glassy phase. From Eq.(H), 
for A > 1 and g > 1/2 we deduce a sort of duality relation, 

Z{u;,Xg) = \Z{u;/VX,g) (13) 

connecting high frequencies, large disorder, to lower frequencies and smaller 
disorder. Z is inserted into the global density-of-states p = po + Pi + P2 
and the indicated integrals performed, with the results shown in Fig. 3. For 
g > 1/2, two fixed points in p are discerned: one just below the characteristic 
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frequency uq and the other just below Suo- A new high-frequency "tail" grows 
beyond the second fixed point, indicating an accumulation of spectral density 
at the highest frequancy as well as at the lowest. 

It is also clear that the two van Hove cusps aXuj/ujQ — -\/2 and 2 are sharp 
only at g — Q-i and that p losses its van Hove cusps and becomes analytic in 
uj for any nonzero value of the disorder, i.e. already in the disordered crystal 
and not just in the glassy phase! The theory works extraordinarily well for g 
from zero into the glassy phase at 0.5 + e, with almost every noncrystalline 
material being excellently characterized at low T by some e < 0.1, and it 
remains trustworthy in every detail until g exceeds 1, when the mean- field 
approximation to the quartic terms fails at high frequencies. 

THE SPECIFIC HEAT. With a knowledge of p the model's ther- 
modynamic properties, including specific heat, entropy and the other such 
functions are readily obtained. For g < 1/2, we find at low T, 



NkB \%J 



where B{g) = (47r2v^/5) ( 1 + 7( [ (1/2) + ^(1/4) -^^ ]-3/2 _ i ) ) and 
7 1.43. For 5- > 1/2, 

c^^-^^Aig) (I-) (15) 



where A{g) = 7r{2/3f/yW[3, 1](1 - {2g)-^/^) ^ 1.404(1 - {2g)-''/^). 

The high-T specific heat is reasonably universal, always tending smoothly 
to the Dulong-Petit limit because J^^ dujpo{uj) = 3. Figure 4 illustrates the 
separate contributions Cq, c-i and C2 from the crystal, disorder and anhar- 
monicity respectively, at g = 0.7, over the interval < T/Tq < 0.6. 

We intend to publish the mathematical details of the various calculations 
elsewhere, but in the meantime the interested reader may request them from 
the corresponding author, M.M. 
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Captions List 



Figure 1: Real part of the self-energy Z{uj) as a function of frequency for 
several disorder concentrations, below and above critical {g = 1/2). 

Figure 2: Imaginary part of the self-energy Z{uj) as a function of frequency 
for several disorder concentrations, below and above critical {g = 1/2). 

Figure 3: Total density of states p(co') = po{uj) + pi{uj) + p2{oj) as a function 
of frequency for a wide range of disorder concentrations, proportional to 
(yf^. Arranged according to the height of the central maximum, we have 
= 0, 0.1, 0.2, ... down to = 1. 

Figure 4: The three distinct contributions to the low-temperature specific 
heat from Eqs. (H), and (|) in the glassy phase {g = 0.7). 
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